#!/usr/bin/env python3
# -*- coding: utf-8 -*-

import matplotlib.pyplot as plt
plt.clf()

VE=13.6E-3 #volume à l equivalence en L
Cbase=0.10 #concentration en quantite de matiere du reactif titrant en mol/L
nAH_i=Cbase*VE #quantite de matiere initiale du reactif titre en mol
V=[i*VE/13.6 for i in range(25)] #crée 25 valeurs de volume de solution titrante ajoute
#valeurs rangees dans une liste
nAH=[]
nOH=[]
nNa=[]
nA=[]

for v in V :
    if v<VE :
        nAH.append(nAH_i-Cbase*v)
        nOH.append(0)
        nNa.append(Cbase*v)
        nA.append(Cbase*v)
    
    else :
        nAH.append(0)
        nOH.append(Cbase*(v-VE))
        nNa.append(Cbase*v)
        nA.append(Cbase*VE)

plt.plot(V, nAH, 'b+-', label="AH")
plt.plot(V, nOH, 'go-', label="OH-")
plt.plot(V, nNa, 'rx', label="Na+")
plt.plot(V, nA, 'cs', label="A-")
plt.xlabel("Volume de solution titrante en L")
plt.ylabel("Quantite de matiere en mol")
plt.legend()
plt.show()